clc,clear
p=[0.8 0.1 0.1;0.5 0.1 0.4;0.5 0.3 0.2];
p=sym(p'); %把p的转置矩阵变成符号矩阵，为了求精确解
[v,d]=eig(p) %求符号矩阵的特征向量矩阵v,特征值矩阵d
d=diag(d); d=double(d); %提取特征值，为了比较大小，把符号值转换成double类型
ind=find(d==max(d));  %求最大特征值的地址
p=v(:,ind)/sum(v(:,ind))  %把最大特征值对应的特征向量化成概率向量
